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Abstract 

Various perturbative and non-perturbative many-body techniques are discussed in this work. 
Especially, we will focus on the summation of so-called Parquet diagrams with emphasis on 
applications to finite nuclei. Here, the subset of two-body Parquet equations will be discussed. A 
practical implementation of the corresponding equations for studies of effective interactions for 
finite nuclei is outlined. 
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1 Introduction 



Traditional shell-model studies have recently received a renewed interest through large 
scale shell model calculations in both the lpOf shell and the 2sldQg-i/2 shells with the 
inclusion of the O/in/2 intruder state as well. It is now therefore fully possible to perform 
large-scale shell-model investigations and study the excitation spectra for systems with 
some 10 million basis states. With recent advances in Monte Carlo methods one is also 
able to enlarge the dimensionality of the systems under study considerably, and important 
information on e.g., ground state properties has thereby been obtained. 
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An important feature of such large scale calculations is that it allows one to probe the 
underlying many-body physics in a hitherto unprecedented way. The crucial starting point 
in all such shell-model calculations is the derivation of an effective interaction, be it either 
an approach based on a microscopic theory starting from the free nucleon-nucleon (NN) 
interaction or a more phenomenologically determined interaction. In shell-model studies 
of e.g., the Sn isotopes, one may have up to 31 valence particles or holes interacting 
via e.g., an effective two-body interaction. The results of such calculations can therefore 
yield, when compared with the availiable body of experimental data, critical inputs to 
the underlying theory of the effective interaction. Until very recently, realistic shell-model 
effective interactions have mainly been applied to nuclei with two or a few valence particles 
beyond closed shells, such as the oxygen and calcium isotopes. Thus, by going to the tin 
isotopes, in which the major neutron shell between neutron numbers 50 and 82 is being 
filled beyond the 100 Sn closed shell core, we have the opportunity of testing the potential of 
large-scale shell-model calculations as well as the reliability of realistic effective interactions 
in systems with many valence particles. 

Clearly, although the NN interaction is of short but finite range, with typical interparticle 
distances of the order of 1 ~ 2 fm, there are indications from both studies of few-body 
systems such as the triton and infinite nuclear matter, that at least three-body interactions, 
both real and effective ones, may be of importance. Thus, with many valence nucleons 
present, such large-scale shell-model calculations may tell us how well e.g., an effective 
interaction which only includes two-body terms does in reproducing properties such as 
excitation spectra and binding energies. 

This work deals therefore with various ways of deriving the effective interaction or effec- 
tive operator needed in shell-model calculations, starting from the free NN interaction. 
Normally, the problem of deriving such effective operators and interactions are solved in 
a limited space, the so-called model space, which is a subspace of the full Hilbert space. 
The effective operator and interaction theory is then introduced in order to systematically 
take into account contributions from the complement (the excluded space) of the chosen 
model space. Several formulations for such expansions of effective operators and interac- 
tions exit in the literature, following time-dependent or time-independent perturbation 
theory Formulations like the coupled-cluster method or exponential ansatz |5|,|7|-|13|, 

the summation of the Parquet class of diagrams |I^-pT|], or the so-called Q-box method 



with the folded-diagram formulation of Kuo and co-workers [||,[|] have been extensively 
applied to systems in nuclear, atomic, molecular and solid-state physics. Here we will focus 
on the above-mentioned Q-box approach and the summation of the so-called Parquet dia- 
grams. For the description of other many-body methods such as the Hypernetted-chain or 
the correlated basis function methods, monte-carlo related methods |2l||2"4fl , Unitary- 
correlation operator method |J etc. see other contributions in this volume. 

The Q-box has been introduced in Rayleigh-Schrodinger perturbation theory as the defi- 
nition of all non-folded diagrams to a given order in the expansion parameter, in nuclear 
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physics the so-called G-matrix. The G-matrix renders the free NN interaction V, which is 
repulsive at small internucleon distances, tractable for a perturbative analysis through the 
summation of ladders diagrams to infinite order. Stated differently, the G-matrix, through 
the solution of the Bethe-Brueckner-Goldstone equation, accounts for the short-range cor- 
relations involving high-lying states. Folded diagrams are a class of diagrams which arise 
due to the removal of the dependence of the exact model-space energy in the Brillouin- 
Wigner perturbation expansion. Through the Q-box formulation and its derivatives, this 
set of diagrams can easily be summed up. 



In addition to the evaluation of folded diagrams and the inclusion of ladder diagrams to 
infinite order included in the G-matrix, there are other classes of diagrams which can be 
summed up. These take into account the effect of long-range correlations involving low- 
energy excitations. A frequently applied formalism is the Tamm-Dancoff (TDA) or the 
random-phase (RPA) approximations. In their traditional formulation one allows for the 
summation of all particle-hole excitations, both forward-going and backward going inser- 
tions. This set of diagrams, as formulated by Kirson |25|] and reviewed in Ref . |26f , should 
account for correlations arising from collective particle-hole correlations. Another possibil- 
ity, is to include any number of particle-particle and hole-hole correlations in the screening 
of particle-hole correlations. The inclusion of these kind of correlations is conventionally 
labelled particle-particle (pp) RPA. It has been used both in nuclear matter studies |27| - 
30| and in evaluations of ground state properties of closed-shell nuclei pT|-|33|j. Recently, 



Ellis, Mavromatis and Miither [|32| , |33|1 have extended the pp RPA to include the particle- 
hole (ph) RPA, though screening of two-particle-one-hole (2plh) and two-hole-one-particle 
(2hlp) vertices was not included. The latter works can be viewed as a step towards the 
full summation of the Parquet class of diagrams. The summation of the Parquet diagrams 
entails a self-consistent summation of both particle-particle and hole-hole ladder diagrams 
and particle-hole diagrams. Practical solutions to this many-body scheme for finite nuclei 
will be discussed here. 



This work falls in six sections. In the next section we present various definitions pertinent to 
the determination of effective interactions, with an emphasis on perturbative methods. The 
resummation of the ladder type of diagrams is then presented in section |||. In that section 
we also discuss the summation of so-called folded diagrams which arise in the evaluation 
of valence space effective interactions. Further perturbative corrections are also discussed 
and selected results for light nuclei in the IsOd and lpO/ shells and for Sn isotopes are 
presented. 

The summation of the TDA and RPA class of diagrams is discussed in section |j. Other 
screening corrections are also discussed in that section. The self-consistent approach to 
the summation of both ladder type diagrams and screening terms through the solution 
of equations for the Parquet class of diagrams will be presented in section [5| Concluding 
remarks are given in section |^. 
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2 Perturbative methods 



In order to derive a microscopic approach to the effective interaction and/or operator 
within the framework of perturbation theory, we need to introduce various notations and 
definitions pertinent to the methods exposed. In this section we review how to calculate 
an effective operator within the framework of degenerate Rayleigh-Schrodinger (RS) per- 
turbation theory P||. 

It is common practice in perturbation theory to reduce the infinitely many degrees of 
freedom of the Hilbert space to those represented by a physically motivated subspace, the 
model space. In such truncations of the Hilbert space, the notions of a projection operator 
P onto the model space and its complement Q are introduced. The projection operators 
defining the model and excluded spaces are defined by 

i=l 

and 

oo 

Q= I^X^I, (2) 

i=D+l 

with D being the dimension of the model space, and PQ = 0, P 2 = p, Q 2 = Q and 
P + Q — I. The wave functions are eigenfunctions of the unperturbed hamiltonian 
Hq = T + U, where T is the kinetic energy and U and appropriately chosen one-body 
potential, that of the harmonic oscillator (h.o.) in most calculations. The full hamiltonian 
is then rewritten as H = H + Hi with Hi = V — U, V being e.g. the NN interaction. The 
eigenvalues and eigenfunctions of the full hamiltonian are denoted by |^ a ) and E a , 

H\* a ) =E a \* a ). (3) 

Rather than solving the full Schrodinger equation above, we define an effective hamiltonian 
acting within the model space such that 

PH eS P |* a ) = E a P |¥ a ) = E a |$ a > (4) 

where |$ a ) = P |\P a ) is the projection of the full wave function onto the model space, the 
model space wave function. In RS perturbation theory, the effective interaction H e $ can 
be written out order by order in the interaction Hi as 

PH cS P = PHiP + PHi^-HiP + PHi^-Hi^-HiP + . . . , (5) 

e e e 

where terms of third and higher order also include the aforementioned folded diagrams. 
Further, e = uj — H , where ui is the so-called starting energy, defined as the unperturbed 
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energy of the interacting particles.. Similarly, the exact wave function \^ a ) can now be 
written in terms of the model space wave function as 

|#a> = |$«) + jHt + jHt^Hx |$ a ) + . . . (6) 
The wave operator is often expressed as 

n = l + X, (7) 

where \ is known as the correlation operator. The correlation operator generates the com- 
ponent of the wave function in the Q-space and must therefore contain at least one pertur- 
bation. Observing that PQP = P, we see that the correlation operator x has the properties 

P X P = 0, QQP = Q X P = XP- (8) 

Since = Q ^f 1 ^ determines the wave operator only when it operates to the right on 
the model space, i.e., only the QP part is defined, the term QQ never appears in the theory, 
and we could therefore add the conditions QxQ = and PxQ = to Eq. (H). This leads 
to the following choice for x 

x = QxP- (9) 

This has been the traditional choice in perturbation theory JTIJ 



The wave operator Q can then be ordered in terms of the number of interactions with the 
perturbation H\ 

n = i + nw + ^ (2) + (io) 

where Vt^ means that we have n Hi terms. Explicitely, the above equation reads 

where e are the unperturbed energies of the P-space and Q-space states defined by H . The 
greek letters refer to Q-space states, whereas a latin letter refers to model-space states. The 
second term in the above equation corresponds to fl^ 1 ' while the third and fourth define f^ 2 ). 
Note that the fourth term diverges in case we have a degenerate or nearly degenerate model 
space. It is actually divergencies like these which are to be removed by the folded diagram 
procedure to be discussed in the next section. Terms like these arise due to the introduction 
of an energy independent perturbative expansion. Conventionally, the various contributions 
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to the perturbative expansion are represented by Feynman-Goldstone diagrams. In Fig. [1] 
we display the topologically distinct contributions to first order in the interaction V to 
the wave operator Eq. ([12]). These diagrams all define the correlation operator x to first 
order in V. Diagrams with Hartree-Fock contributions and single-particle potential terms 
are not included. The possible renormalizations of these diagrams will be discussed in the 
next three sections. The reader should note that with respect to the nomenclature in Eq. 
(|T2|), we will hereafter employ the following notation in our discussion of various diagrams 
and vertex renormalizations: 

- roman letters p, q, r,s,t, . . . refer to particle single-particle states, either within the 
model-space or from the excluded space. An arrow pointing upwards represents such 
a particle state. 

- greek letters a, (5, 7, 8 . . . refer to hole single-particle states. An arrow pointing downward 
is a hole state. 



(a) 





5 






/ p\ 






(b) 






r 


y 


j 





(c) 



(d) 



(e) 



(f) 



(9) 



Fig. 1. The various vertices to first order in the interaction V which contribute to the wave 
operator £1 = 1 + x- Hartree-Fock terms are not included. Possible hermitian conjugate diagrams 
are also not shown. Indicated are also possible angular momentun coupling orders. 



2. 1 Expressions for the wave operator 



We end this section with the equations for the diagrams in Fig. [I] representing x to first 
order in V. Moreover, in order to introduce the various channels needed to sum the Parquet 
class of diagrams, we will find it convenient here to classify these channels in terms of 
angular momentum recouplings. Later on, we will also introduce the pertinent definitions 
of energy and momentum variables in the various channels. The nomenclature we will follow 
in our labelling is that of Blaizot and Ripka, see Ref. |18| chapter 15. All matrix elements 
in the definitions below are antisymmetrized and unnormalized. The first channel is the 
[12] channel, or the s-channel in field theory, and its angular momentum coupling order 
is depicted in Fig. 0. In this figure we do not distinguish between particles and holes, all 
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J 



(a) (b) (c) 

Fig. 2. Coupling order for the [12] (a), [13] (b) and [14] (c) channels. 

single-particle labels 1, 2, 3, 4 can represent either a hole or particle single-particle state. It 
is the coupling order which is of interest here. The matrix element V ^ in the [12] channel 
is then 

V^L=<(12)J|y|(34)J), (12) 

meaning that the single-particle state 1(3) couples to the state 2(4) to yield a total angular 
momentum J. This channel is commonly denoted as the particle-particle (pp) or particle- 
particle-hole-hole (pphh) channel, meaning that when we will sum classes of diagrams to 
infinite order in this channel, the only intermediate states which are allowed are those of 
a pphh character, coupled to a final J in the above order. In the next section we will 
explicitely discuss ways to evaluate the equations for the [12] channel. This coupling order 
is also the standard way of representing effective interactions for shell-model calculations. 
If we now specialize to particles and holes (these matrix elements were shown in Fig. p]) 
we obtain for the case with particle states only, diagram (a), 

V 2p = vWj=((pq)J\V\(rs)J). (13) 
The corresponding one for holes only, diagran (b), is 

V* = v!$ 6J =(( a l3)J\V\bS)J). (14) 

Thus, in the forthcoming discussion, we will always employ as our basic notation for a 
matrix element that of the [12] channel, meaning that matrix elements of the other two 
channels can always be rewritten in terms of those in [12] channel We see this immediately 
by looking at the expression for the matrix element in the [13] channel, the t-channel in 
field theory, see Fig. 0(b). Here the single-particle state 3(4) couples to the single-particle 



state 1(2)111. Through simple angular momentum algebra we have 
J' 32 k J' 



1 In a Goldstone-Feynman diagram in an angular momentum representation, the coupling di- 
rection will always be from incoming single-particle states to outgoing single-particle states. 



S 



where the symbol with curly brackets represents a 6j-symbol and J' = \ / 2J' + 1. In a 
similar way we can also express the matrix element in the [14] channel, the w-channel in 
field theory, through 



VS&J = U~) h+ji+J+2h J' 2 { H 31 3 \ (16) 

h ja J' 

It is also possible to have the inverse relations or to express e.g., the [14] channel through 
the [13] channel as 

vSSlj = U-) 2h+2h+2j3 J' 2 l 34 31 ' \ V&ij. (17) 
J' h h -J' 



The matrix elements defined in Eqs. (p^)-([T7|) and the inverse relations are the starting 
points for various resummation of diagrams. In the next section we will detail ways of 
solving equations in the [12] channel, whereas various approximations for the [13] channel 
and [14] channel such as the TDA and RPA and vertex and propagator renormalization 
schemes will be discussed in section |j. Finally, how to merge self-consistently all three 
channels will be discussed in section ^. 

We end this section by giving the expressions in an angular momentum basis i for the 
remaining diagrams of Fig. [I]. The coupling order is indicated in the same figure. 

Thus, the 2plh vertex V^pih, diagram (c) in Fig. [I], is coupled following the prescription of 
the [13] channel and reads 



> / 2plh — V pqr a J ~ Z^V ) \ I WaJ'- 

J' ( jq ja J' J 



The 2p2h ground-state correlation V^h, diagram (d), which will enter in the RPA sum- 
mation discussed in section [| is given by, the coupling order is that of the [13] channel, 



^ P2h = O = j2(-y^ +J+J 'j' 2 \ 3a3pJ \ v ™ Jt . do) 

J' jq J (3 J' 



2 Note that we only include angular momentum factors, other factors coming from the diagram 
rulespl], like number of hole lines, number of closed loops etc. are omitted here. 
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The 2hlp vertex V^hip, diagram (e), still in the representation of the [13] channel, is defined 

as 



J' I J/3 J P J' J 



= ^ = E(-)" +i ' +J+J '> Jl M V2&,. (20) 



Note well that the vertices of Eqs. (]T8|)-(|20|) and their respective hermitian conjugates 
can all be expressed in the [14] channel or [12] channel as well. However, it is important 
to note that the expressions in the various channels are different, and when solving the 
equations for the various channels, the renormalizations will be different. As an example, 
consider the two particle- hole vertices Vph of Fig. [j], i.e., diagrams (f) and (g). Diagram 
(g) is just the exchange diagram of (f) when seen in the [12] channel. However, if (f) is 
coupled as in the [13] channel, recoupling this diagram to the [14] channel will not give two 
particle-hole two-body states coupled to a final J but rather a particle-particle two-body 
state and a hole-hole two-body state. But why bother at all about such petty details? The 
problem arises when we are to sum diagrams in the [13] channel and [14] channel. In the 
[12] channel we allow only particle-particle and hole-hole intermediate states, whereas in 
the [13] channel and [14] channel we allow only particle-hole intermediate states, else we 
may risk to double-count various contributions. If we therefore recouple diagram (f ) to the 
[14] representation, this contribution does not yield an intermediate particle-hole state in 
the [14] channel. Thus, diagram (f), whose expression is 



v = vsl = u-y^ +j+j, j' \ jajp \\ vz aJ „ (2D 



u-y^ +j+j, j' 2 l iaip J )v£?u 

J ' ( 3P 3q ^ J 



yields a particle-hole contribution only in the [13] channel, whereas the exchange diagram 
(g), which reads 

Vp- = vXIj = £(-)*«-■»**' ( 3 " * J \ V*JU (22) 

J' { 30 3q J' 

results in the corresponding particle-hole contribution in the [14] channel. In electron gas 
theory, the latter expression is often identified as the starting point for the self-screening 
of the exchange term. In the discussion of the TDA series in section |4] we will give the 
expressions for the screening corrections based on Eqs. (|21|) and (p2|) . 



An important aspect to notice in connection with the latter equations and the discussions 
in sections ^ and [5] is that 

y M _ _y[13] (no) 
V pf3qaJ — v pf3aqJ' l zc V 
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i.e., just the exchange diagram, as it should be. This is however important to keep in mind, 
since we later on will sum explicitely sets of diagrams in the [13] channel and the [14] 
channel, implying thereby that we will obtain screening and vertex corrections for both 
direct and exchange diagrams. 



3 Summation of diagrams in the [12] channel 



In order to write down the equation for the renormalized interaction in the [12] channel 
we need first to present some further definitions. We will also assume that the reader has 
some familiarity with the theory of Green's function. Thorough discussions of such topics 
can be found in the recent reviews of Dickoff and Miither |55 and Kuo and Tzeng p6| . 



In our presentation below we will borrow from these works and the monograph of Blaizot 
and Ripka |18|. The vertex r' 12 ' is in lowest order identical with the interaction and 



obeys also the same symmetry relations as V, i.e., 

r [i2] _ r [i2] _ _ F [i2] _ p[12] m x 

1 1234J — 1 2143J ~~ 1 2134J — 1 1243J- \^) 

We also need to define energy variables. Since we are going to replace the interaction V 
with the G-matrix, or certain approximations to it, defined below in all of our practical 
calculations, the momentum variables are already accounted for in G. The basis will be 
that of harmonic oscillator wave functions, and the labels 1234 will hence refer to oscillator 
quantum numbers, which in turn can be related to the momentum variables. The labels 
1234, in addition to representing single-particle quantum numbers, define also the energy 
of the single-particle states. With a harmonic oscillator basis, the starting point for the 
single-particle energies £1,2,3,4 are the unperturbed oscillator energies. When iterating the 
equations for I^ 12 ', self-consistent single-particle energies can be introduced. The total 
energy in the [12] channel s is 

s = Si + e 2 = e 3 + £4. (25) 

The equation for the vertex r' 12 ] is, in a compact matrix notation, given by [1^ 



r [12] = y[12] + y [12] Wr [12] j (26) 

where g is the one-body Green's function representing the intermediate states. The dia- 
grammatic expression for this equation is given in Fig. |3|. The expression of Eq. ( p6|) is 
known as the Feynman-Galitskii equation. This equation is normally solved iteratively. In 
the first iteration the irreducible vertex yt 12 ! is then often chosen as the bare NN inter- 
action. This interaction is then typically assumed to be energy independent and we can 
drop the s dependence of V^ 12 ]. Moreover, the single-particle propagators are chosen as the 
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V 






r 








Fig. 3. (a) represents the two-body vertex T function while (b) represents the self-energy E. 
unperturbed ones. The first iteration of Eq. (^6j) can then be rewritten as 



r [i2] r x y[i2] ,ivi/ [121 c' 12 'r [121 (s) 

1 1234jl 6 J _ > / 1234J + 9 L v 1256jy 1 5634jl 6 J 
Z 56 



(27) 



with the unperturbed particle-particle and hole-hole propagators 



;[i2] 



Vhh 



s — e 5 — e 6 + ir] s - e 5 - e 6 - vr] ' 



(28) 



which results from the integration over the energy variable in the product of the two 
single-particle propagators in Eq. (|2"rj|). The factor 1/2 follows from one of the standard 
Goldstone-Feynman diagram rules f34|, which state that a factor 1/2 should be associated 
with each pair of lines which starts at the same interaction vertex and ends at the same 
interaction vertex. The reader should note that the intermediate states 56 can represent 
a two-particle state or a two-hole state. In Eq. (|28|) we have assumed unperturbed single- 
particle energies. In our iterations we will approximate the single-particle energies with their 
real part only. Thus, the two-particle propagator with renormalized single-particle 
energies has the same form as the unperturbed one. The operators Qp^ and Qy£ ensure 
that the intermediate states are of two-particle or two-hole character. In order to obtain 
a self-consistent scheme, Eq. fl2"T| ) has also to be accompanied with the equation for the 
single-particle propagators g given by Dyson's equation 



9 = go - go^g, 



(29) 



with g being the unperturbed single-particle propagator and £ the self-energy. We will 
however defer a discussion of these quantities to section [|. Here it will suffice to state that 
the self-energy is related to the vertex as 



(30) 
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The similarity sign is meant to indicate that, although being formally correct, great care 
has to be exercised in order not to double-count contributions to the self-energy [[L7 |. The 
set of equations for the vertex function and the self-energy is shown pictorially in Fig. [| 
Assume now that we have performed the first iteration. The question which now arises 
is whether the obtained vertex r' 1 ^ from the solution of Eq. fl2"T|) should replace the bare 
vertex in the next iteration. Before answering this question, let us give some examples 
of diagrams which can be generated from the first iteration. These contributions are given 
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(c) 



(d) 



(e) 



(f) 



Fig. 4. Diagrams (a)-(d) give examples of diagrams which are summed up by the use of Eq. (P6|). 
Diagrams (e) and (f) are examples of core-polarization terms which are not generated by the [12] 
channel. 

by e.g., diagrams (a)-(d) in Fig. |j. Diagrams (a) and (b) are examples of contributions to 
second order in perturbation theory, while diagrams (c) and (d) are higher order terms. 
Diagrams (e) and (f) are higher-order core-polarization terms, which can e.g., be generated 
through the solution of the equations for the [13] and [14] channels discussed in the next 
section. If we were to include diagrams (a)-(d) in the definition of the bare vertex in our next 
iteration, we are prone to double-count, since such contributions are generated once again. 
Diagrams which contain particle-hole intermediate state are however not generated by the 
solution of Eq. (|27|). We need therefore to define the vertex V^ 12 l used in every iteration 
to be the sum of diagrams which are irreducible in the [12] channel. With irreducible we 
will mean all diagrams which cannot be reduced to a piece containing the particle states 
12 entering or leaving the same interaction vertex and another part containing the states 
34 at the same interaction vertex by cutting two internal lines. Clearly, if we cut diagrams 
(a) and (b) we are just left with two bare interaction vertices. Similarly, cutting two lines 
of an intermediate state in diagrams (c) and (d) leaves us with two second-order terms of 
the type (a) and (b) and one bare interaction. Diagrams (e) and (f) are however examples 
of diagrams which are irreducible in the [12] channel. Diagram (e) is irreducible in the [13] 
channel, but not in the [14] channel. Similarly, diagram (g) is reducible in the [13] channel 
and irreducible in the [14] channel. This means that, unless we solve equations similar to 



13 



Eq. ( p7|) in the [13] channel and [14] channels as well, changes from further iterations of 



In the remaining part of this section, we will try to delineate ways of solving the above 
equations, and discuss possible approximations, their merits and faults. First of all, we will 
reduce the propagator of Eq. (|28| ) to only include particle-particle intermediate states. This 
will lead us to the familiar G-matrix in nuclear many-body theory. Based on the G-matrix, 
we will construct effective interactions through perturbative summations. Applications of 
such effective interacions to selected nuclei will then be discussed. Thereafter, we will try to 
account for hole-hole contributions and end this section with a discussion on self-consistent 
determinations of the single-particle energies through the solution of Dyson's equation. 



3.1 The G -matrix 



In nuclear structure and nuclear matter calculations one has to face the problem that any 
realistic NN interaction V exhibits a strong short-range repulsion, which in turn makes a 
perturbative treatment of the nuclear many-body problem prohibitive. If the interaction has 
a so-called hard core, the matrix elements of such an interaction V evaluated for an 
uncorrelated two-body wave function ip(r) diverge, since the uncorrelated wave function is 
different from zero also for relative distances r smaller than the hard-core radius. Similarly, 
even if one uses interactions with softer cores, the matrix elements of the interaction become 
very large at short distances. The above problem was however overcome by introducing the 
reaction matrix G^ (displayed by the summation of ladder type of diagrams in Fig. |] with 
particle-particle intermediate states only), accounting thereby for short-range two-nucleon 
correlations. The G^-matrix represents just a subset to the solution of the equations for 
the interaction r' 12 ] in the [12] channel, we have clearly neglected the possibility of having 
intermediate states which are of the hole-hole type. The matrix elements of the interaction 
then become 



where \1/ is now the correlated wave function containing selected correlation from the ex- 
cluded space. By accounting for these correlations in the two-body wave functon the ma- 
trix elements of the interaction become finite, even for a hard-core interaction V. Moreover, 
as will be discussed below, compared with the uncorrelated wave function, the correlated 
wave function enhances the matrix elements of V at distances for which the interaction is 
attractive. The type of correlations which typically are included in the evaluation of the 
G^'-matrix are those of the two-particle type. If we label the operator Q in this case by 





(31) 
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Qpp i we can wr it e the integral equation for the G- matrix as 



G [12U s) = V[ 12] + V 112] ^ G [12]/ s) (32) 

s- Ho + tT) K " 

implicitely assuming that limi] — > oo. The variable s represents normally the unperturbed 
energy of the incoming two-particle state. We will suppress vq in the following equations. 
Moreover, since one is often interested only in the Gt 12 l-matrix for negative starting energies, 
the Gt 12 ]-matrix commonly used in studies of effective interactions has no divergencies. 

We can also write 

G™(s) = V™ + V^Q^ 1 QMGM(s). (33) 

s — Qpp HqQpp 

The former equation applies if the Pauli operator commutes with the unperturbed 
hamiltonian Hq, whereas the latter is needed if [Hq, Q^} ^ 0. Similarly, the correlated 
wave function \l/ is given as 

m = 1^ + ^-0^^), (34) 



Ho 



or 



s — Qpp ri ovpp 

In order to evaluate the G' 12 '-matrix for finite nuclei, we define first a useful identity 
following Bethe, Brandow and Petschek ||37|| . Suppose we have two different G-matricesH, 
defined by 

Gl = Vl + Vl 9l Gu ( 36 ) 
ex 

and 

G 2 = V 2 + V 2 —G 2 , (37) 

e 2 

where Q\je\ and Q 2 /e 2 are the propagators of either Eq. fl32|) or Eq. (|33|) . Gi and G 2 are 
two different G-matrices having two different interactions and/or different propagators. 



3 For notational economy, we drop the superscript I 12 !. Furthermore, in the subsequent discussion 
in this subsection it is understood that all operators Q refer to particle-particle intermediate states 
only. The subscript pp is also dropped. 
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We aim at an identity which will enable us to calculate G\ in terms of G2, or vice versa. 
Defining the wave operators 

n 1 = l + Q± Gl , (38) 
ei 

and 

n 2 = 1 + ^G 2 , (39) 
e 2 

we can rewrite the above G-matrices as 

Gi = VtQx, (40) 

and 

G 2 = v 2 n 2 . (41) 
Using these relations, we rewrite G\ as 



g 1 =g 1 - g\ (n x -i- ^d) + (n\ - 1 - g 2 ^) d 

Gj + G \ 2 (9l _ ^ G + nl Gl _ Gt n (42) 
V ei e 2 / 

and using Eqs. ( ^0|) and (|4l| ) we obtain the identity 

Gi = G\ + f 5l -9l) Gl + nj(y _ % ) ni . (43) 
V ei e 2 / 

The second term on the rhs. is called the propagator-correction term; it vanishes if G\ 
and G 2 have the same propagators. The third term is often referred to as the potential- 
correction term, and it disappears if G\ and G 2 have the same potentials. The reader may 
now ask what is the advantage of the above identity. If we assume that by some physical 
reasoning we are able to calculate G 2 and that the expression for G 2 can be calculated 
easily, and further that G 2 is a good approximation to the original G-matrix, then we can 
use the above identity to perform a perturbative calculation of G\ in terms of G 2 . 

Before we proceed in detailing the calculation of the G-matrix appropriate for finite nuclei, 
certain approximations need be explained. 

As discussed above, the philosophy behind perturbation theory is to reduce the intractable 
full Hilbert space problem to one which can be solved within a physically motivated model 
space, defined by the operator P. The excluded degrees of freedom are represented by 
the projection operator Q. The definition of these operators is connected with the nuclear 
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Orbit number 



4s-3d-2g-li-0k 


37-45 


3p-2f-lh-0j 


29-36 


3s-2d-lg-0i 


22-28 


2p-lf-0h 


16-21 


2s-ld-0g 


11-15 


lp-Of 


7-10 


ls-Od 


4-6 


Op 


2-3 


Os 


1 



Fig. 5. Classification of harmonic oscillator single-particle orbits. 

system and the perturbative expansions discussed in section |2|. Consider the evaluation 
of the effective interaction needed in calculations of the low-lying states of 18 0. From 
experimental data and theoretical calculations the belief is that several properties of this 
nucleus can be described by a model space consisting of a closed 16 core (consisting of 
the filled Os- and Op-shells) and two valence neutrons in the lsOrf-shell. In Fig. [5] we exhibit 
this division in terms of h.o. sp orbits. The active sp states in the IsOrf-shell are then 
given by the 0d 5 / 2 , 0d 3 / 2 and lsi/2 orbits, labels 4 — 6 in Fig. |5j. The remaining states 
enter the definition of Q. Once we have defined P and Q we proceed in constructing the 
G-matrix and the corresponding perturbative expansion in terms of the G-matrix. There 
are however several ways of choosing Q. A common procedure is to specify the boundaries 
of Q by three numbers, ni, n 2 and n 3 , explained in Fig. |[ For 18 we would choose 



n3 



n2 



0=1 



nl 



Q=0 



t t 

nl n2 



q3 



Fig. 6. Definition of the P (shaded area) and Q operators appropriate for the definition of the 
G-matrix and the effective interaction. Outside the shaded area limited by the boundaries n\, n 2 
and n-3 P = and Q = 1. 



(ni 



3, n 2 = 6,n 3 = oo). Our choice of P-space implies that the single-particle states 
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outside the model space start from the lpO/-shell (numbers 7-10 in Fig. |5|), and orbits 
1, 2 and 3 are hole states. Stated differently, this means that Q is constructed so as to 
prevent scattering into intermediate two-particle states with one particle in the 0s- or 
Op-shells or both particles in the IsOrf-shell. This definition of the Q-space influences the 
determination of the effective shell-model interaction. Consider the diagrams displayed in 
Fig. [7]. Diagram (a) of this figure is just the G-matrix and is allowed in the definition of 













pi' 












pi' 



(a) 



(b) 




AAAA 



Vaaa/ 



(d) 



Fig. 7. Examples of diagrams which may define the effective valence space interaction. The wavy 
line is the G-matrix. 

the Q-box. With our choice (n\ = 3,n 2 = 6,n 3 = oo), diagram (b) is not allowed since 
the intermediate state consists of passive particle states and is already included in the 
evaluation of the G-matrix. Similarly, diagram (c) is also not allowed whereas diagram 
(d) is allowed. Now an important subtlety arises. If we evaluate the G-matrix with the 
boundaries {n\ = 3,n 2 = 10, n 3 = oo), and define the P-space of the effective interaction 
by including orbits 4 to 6 only, then diagrams (b) and (c) are allowed if 7 < pi,p 2 < 10 
In this way we allow for intermediate two-particle states as well with orbits outside the 
model-space of the effective interaction. The reader should notice the above differences, 
i.e. that the Q-space defining the G-matrix and H g r may differ. In order to calculate the 
G-matrix we will henceforth employ a so-called double-partitioned scheme. Let us be more 
specific and detail this double-partitioned procedure. We define first a reference G-matrix 
G in terms of plane wave intermediate states only, meaning that Hq is replaced by the 
kinetic energy operator T only while G has harmonic oscillator intermediate states (this 
is one possible choice for U). We divide the exclusion operator into two parts, one which 
represents the low- lying states Qi and one which accounts for high- lying states Qh, viz. 

Q = Qi + Qh = Qi + Q- 

If we consider 18 O as our pilot nucleus, we may define Qi to consist of the sp orbits of 
the lpO /-shell, orbits 7 — 10 in Fig. [5|, described by h.o. states. Qh represents then the 
remaining orthogonalized intermediate states. Using the identity of Bethe, Brandow and 
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Petschek [|37| of Eq. (|43|) we can first set up G as 



G = V + V-^-G, (44) 
s — 1 

and express G in terms of G as 

G = G + G f — %-] G, (45) 

and we have assumed that G is hermitian and that [Qi,H Q ] = 0. Thus, we first calculate 
a "reference" G- matrix (G in our case), and then insert this in the expression for the full 
G-matrix. The novelty here is that we are able to calculate G exactly through operator 
relations to be discussed below. In passing we note that G depends significantly on the 
choice of Hq, though the low-lying intermediate states are believed to be fairly well repre- 



sented by h.o. states. Also, the authors of ref. [38 demonstrate that low-lying intermediate 



states are not so important in G-matrix calculations, being consistent with the short-range 
nature of the NN interaction. Since we let Qi to be defined by the orbits of the lpO /-shell, 
and the energy difference between two particles in the s<i-shell and pf shell is of the order 
— 14 MeV, we can treat G as a perturbation expansion in G. Eq. fl45[) can then be written 

as 

G = G + g(-%-)g + g(-%-)g (-%-)&+... (46) 

The only intermediate states are those defined by the lpO /-shell. The second term on the 
rhs. is nothing but the second-order particle-particle ladder. The third term is then the 
third-order ladder diagram in terms of G. As shown by the authors of ref. the inclusion 
of the second-order particle-particle diagram in the evaluation of the Q-box, represents a 
good approximation. The unsettled problem is however how to define the boundary between 
Qi and Q h . 

Now we will discuss how to compute G. One can solve the equation for the G-matrix for 
finite nuclei by employing a formally exact technique for handling Q discussed in e.g., Ref. 
||38|| . Using the matrix identity 

q_1_q = I _ ip . 1 _ pL (47) 
QAQ A A PA-^P A' v ; 

with A = s — T, to rewrite Eq. fljip as@ 

G = G F + AG, (48) 



4 We will omit the label G for the G-matrix for finite nuclei, however it is understood that the 
G-matrix for finite nuclei is calculated according to Eq. ( |44| ) This means that we have to include 
the particle-particle ladder diagrams in the Q-box. 
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where G F is the free G- matrix defined as 

G F = V + V—=G F . (49) 
s — 1 

The term AG is a correction term defined entirely within the model space P and given by 

AG = -v\p ~ 1 ~ P\v. (50) 

Employing the definition for the free G-matrix of Eq. (Tj9f ), one can rewrite the latter 
equation as 

AG = -G F -P^ ^P-G F , (51) 

e P^- 1 + e- 1 G F e- 1 )P e K J 

with e = s — T . We see then that the G-matrix for finite nuclei is expressed as the sum 
of two terms; the first term is the free G-matrix with no Pauli corrections included, while 
the second term accounts for medium modifications due to the Pauli principle. The second 
term can easily be obtained by some simple matrix operations involving the model-space 
matrix P only. However, the second term is a function of the variable 713. The convergence 
in terms of was discussed ad extenso in Ref. Q and we refer the reader to that work. 
The equation for the free matrix G F is solved in momentum space in the relative and centre 
of mass system and thereafter transformed to the relevant expression in terms of harmonic 
ocillator single-particle wavefunctions in the laboratory system. This yields final matrix 
elements of the type 

({ab)J\G\{cd)J) (52) 

where G is the given by the sum G = G F + AG. The label a represents here all the single 
particle quantum numbers n a l a j a . 



3.2 Folded diagrams and the effective valence space interaction 



Here we discuss further classes of diagrams which can be included in the evaluation of 
effective interactions for the shell model. Especially, we will focus on the summations of 
so-called folded diagrams. 

One way of obtaining the wave operator Q is through the generalized Bloch equation given 
by Lindgren and Morrison || 

[Q, H }P = QH&P - xPHittP, (53) 
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which offers a suitable way of generating the RS perturbation expansion. Writing Eq. ([53]) 
in terms of we have 

[Q (1 \H ]P = QH 1 P, (54) 

[Q {2) , H ]P = QHxOPP - Q {1) PHiP, (55) 

and so forth, which can be generalized to 

n-l 

[f2 (n) , H ]P = QHM n ~ 1] P - n {n - m) PHM m ~ l) P. (56) 

m=l 

The effective interaction to a given order can then be obtained from Q( n \ see ||. Another 
possibility is obvioulsy the coupled-cluster method discussed elewhere in this volume. 

Here we will assume that we can start with a given approximation to Q, and through an 
iterative scheme generate higher order terms. Such schemes will in general differ from the 
order-by-order scheme of Eq. ([53]). Two such iterative schemes were derived by Lee and 
Suzuki |39[. We will focus on the folded diagram method of Kuo and co-workers H]. 



Having defined the wave operator Q = 1 + \ (note that Q 1 = 1 — x) with \ given by Eq. 
(0) we can obtain 

QHP - X HP + QH X - XH X = 0- (57) 

This is the basic equation to which a solution to x is to be sought. If we choose to work 
with a degenerate model space we define 

PH P = sP, 

where s is the unperturbed model space eigenvalue (or starting energy) in the degenerate 
case, such that Eq. (^) reads in a slightly modified form (H = Hq + Hi) 

(s - QH Q - QH X Q) X = QH X P - xP^P - xP^Qx, 

which yields the following equation for x 

^ = 7^HQ QH ^-7^HQ x ^ P + PH ^ xP) - (58) 

Observing that the P-space effective hamiltonian is given as 

H cS = PHP + PH X = PH P + V eS ( X ), 
with Kff(x) = PHiP + PH&xP, Eq. @ becomes 

x = 7^HQ QH ^-7^HQ xV ^ < 59) 
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Now we find it convenient to introduce the so-called Q-box, denned as 



Q(s) = PH X P + PH X Q- 



QHQ 



QHiP. 



(60) 



The Q-box is made up of non-folded diagrams which are irreducible and valence linked. 
A diagram is said to be irreducible if between each pair of vertices there is at least one 
hole state or a particle state outside the model space. In a valence-linked diagram the 
interactions are linked (via fermion lines) to at least one valence line. Note that a valence- 
linked diagram can be either connected (consisting of a single piece) or disconnected. In the 
final expansion including folded diagrams as well, the disconnected diagrams are found to 
cancel out || . This corresponds to the cancellation of unlinked diagrams of the Goldstone 
expansion. We illustrate these definitions by the diagrams shown in Fig. §. Diagram (a) 










(a) 



(b) 



(c) 



Fig. 8. Different types of valence-linked diagrams. Diagram (a) is irreducible and connected, (b) 
is reducible, while (c) is irreducible and disconnected. 

is irreducible, valence linked and connected, while (b) is reducible since the intermediate 
particle states belong to the model space. Diagram (c) is irreducible, valence linked and 
disconnected. It is worth noting that general form of the Q-box is the same as that of the 
G-matrix, or the equations of the [12] channel or those of the [13] and [14] channels to be 
discussed in section £2[ In Ref. [§J, the Q-box was defined to be the sum all diagrams to 
third order in the G-matrix. 



Multiplying both sides of Eq. (|5^) with PHi and adding PH\P to both sides we get 

PH X P + PH lX = PH X P + PH X Q ^—QH.P - PH X L— X V^( X ), 

s — QHQ s—QHQ 

which gives 



V eS { X ) = Q{s)-PH 1 



1 



s-QHQ 



(61) 



There are several ways to solve Eq. (|6l|) . The idea is to set up an iteration scheme where we 
determine X n and thus Ke(Xn) from Xn-i and V e s( Xn _i). For the mere sake of simplicity 
we write v}g' = V cS ( Xn ). 
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Let us write Eq. (pip as 

= Q(s) - PH, V — Xn vt 1 



1 

s-~QHQ' 



The solution to this equation can be shown to be [35 

Observe also that the effective interaction is V}^ is evaluated at a given model space 
energy s. If v}^ = 1 , the iteration is said to converge. In the limiting case n — > oo, 
the solution agrees with the formal solution of Brandow [§] and Des Cloizeaux j|0| 



Alternatively, we can generate the contribution from n folds the following way. In an re- 
folded Q-box there are of course n + 1 Q-boxes. The general expression for an n-folded 
Q-box is then 

Q-QjQ + QjQjQ-... = ^——P——P . . . —^PQ, (64) 
where we have the constraints 

777.1 + m 2 + • • • m n — n , 

mi > 1, 

m 2 ,m 3 , ...m n >0, 

and 

rrik < n — k + 1. 

The last restriction follows from the fact that there are only n — k + 1 Q-boxes to the right 
of k th Q-box. Thus, it can at most be differentiated n — k + 1 times. We have inserted 
the model-space projection operator in the above expression, in order to emphasize that 
folded diagrams have as intermediate states between successive Q-boxes only model-space 
states. Therefore, the sum in Eq. (0) includes a sum over all model-space states with the 
same quantum numbers such as isospin and total angular momentum. It is understood 
that the Q-box and its derivatives are evaluated at the same starting energy, which should 
correspond to the unperturbed energy of the model-space state. It is then straightforward 
to recast Eq. into the form of Eq. (|62j). 
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Note that although Q and its derivatives contain disconnected diagrams, such diagrams 
cancel exactly in each order [Q], thus yielding a fully connected expansion in Eq. (|6^). How- 
ever, in order to achieve this, disconnected diagrams have to be included in the definition 
of the Q-box. An example is given by diagram (c) Fig. |8|. Such a diagram will generate a 

contribution to the first fold ^jjQ = —Q j Q which cancels exactly diagram (c) when all 
time-ordered contributions to this diagram are accounted for, see Ref. || for more details. 
It is moreover important to note in connection with the above expansion, that a term like 
F\ = QiQ actually means PQiPQP since the Q-box is defined in the model space only. 
Here we have defined Q m = ^j^-Jr- Due to this structure, only so-called folded diagrams 
contain P-space intermediate states. 

The folded diagram expansion discussed above yields however a non-hermitian effective 
interaction. This may happen even at the level of the G-matrix, or any of the effective 
interactions we will derive in this work, from Parquet theory to the order-by-order pertur- 
bation expansion. 

We will therefore end this section with a way to cure this non-hermiticity. A hermitian 
effective interaction has recently been derived by Ellis, Kuo, Suzuki and co-workers fliyBpif 
through the following steps i. To obtain a hermitian effective interaction, let us define a 
model-space eigenstate \b\) with eigenvalue A as 

\h) = E \4> a ) (65) 

a=l 

and the biorthogonal wave function 

h) = J2t ) \^), (66) 

a=l 

such that 

[hK) = 5 x ». (67) 

The model-space eigenvalue problem can be written in terms of the above non-hermitian 
effective interaction unperturbed wave functions 

E b< i X) (^\H + V + VQ X |V 7 > = E x bi x l (68) 

7=1 

The exact wave function expressed in terms of the correlation operator is 

|* a ) = (1+x)|Va>. (69) 



The reader who wishes more details can consult Refs. |l|,[4T||. 
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The part x \^x) can be expressed in terms of the time-development operator or using the 
time-independent formalism as 

^ = E^QHQ QVPm ' (?0) 

where Q is the exclusion operator. Note that this equation is given in terms of the Brillouin- 
Wigner perturbation expansion, since we have the exact energy E\ in the denominator. 

Using the normalization condition for the true wave function we obtain 

(* 7 |* A > N X 6 X , = (V 7 | (1 + X ] X |-0a> , (71) 

where we have used the fact that (V^l X I^a) = 0. Recalling that the time-development 
operator is hermitian we have that x^X is a l so hermitian. We can then define an orthogonal 
basis d whose eigenvalue relation is 

E<^lx f xl^K = /4< ( 72 ) 



with eigenvalues greater than 0. Using the definition in Eq. (|70|), we note that the diagonal 
element of 

<V>a| X ] X |-0a> = (V>aI pvq 1 QVP |Va) , (73) 

which is nothing but the derivative of the Q-box, with an additional minus sign. Thus, 
noting that if 7 7^ A 

(* 7 |* A > = = (^|Va) + (V 7 I X ] X |-0a> , (74) 
we can write x X m operator form as 

X j X = ~Y, \^o) M Qi(E a ) |^«> {^ a \ - E |^«) (V»o|^> \^p) • (75) 



Using the new basis in Eq. ([72]) , we see that Eq. ( ffl|) allows us to define another orthogonal 
basis h 



/4 + 1 E ^T/flf = E dfip^x, (76) 



where we have used the orthogonality properties of the vectors involved. The vector h was 
used by the authors of Ref. to obtain a hermitian effective interaction as 

W*\ Ks m) = r~ 9 = r^ = > (77) 

lul + 1 + J 4 + 1 
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where (her) and (nher) stand for hermitian and non-hermit ian respectively. This equation 
is rather simple to compute, since we can use the folded-diagram method to obtain the 
non- hermitian part. To obtain the total hermitian effective interaction, we have to add the 
H Q term. The above equation is manifestly hermitian. Other discussion of the hermiticity 
problem can be found in Refs. ]5],|7| . The remaining question is how to evaluate the Q-box. 
Obviously, we are not in the position where we can evaluate it exactly, i.e., to include 
all possible many-body terms. Rather, we have to truncate somewhere. Several possible 
approaches exist, but all have in common that there is no clear way which tells us where to 
stop. However, as argued by the authors of Ref. |T7[] , there is a minimal class of diagrams 
which need to be included in order to fulfull necessary conditions. This class of diagrams 
includes both diagrams which account for short-range correlations such as the G-matrix 
and long-range correlations such as those accounted for by various core-polarization terms. 
The importance of such diagrams has been extensively documented in the literature and 
examples can be found in Refs. PIJZEfl. In Ref. [[| we included all core-polarization con- 
tributions to third-order in the G-matrix, in addition to including other diagrams which 
account for short-range correlations as well. 



In the next subsection we present results from recent large scale shell-model Monte Carlo 
calculations for light nuclei in the IsOd and lpOf shells based on an effective interac- 
tion for the two shells. The effective interaction was derived following the methods outlined 
hitherto in this section, taking into account the non-hermiticity which arises when dealing 
with more than one major shell. The closed shell core is 16 O and all diagrams through 
second-order where employed in the definition of the Q-box. Folded diagrams were cal- 
culated using the abovementioned folded-diagrams method and the interaction was made 
explicitely hermitian through Eq. ([H]). A discussion of results for Sn isotopes will also be 
given. 



3.3 Selected applications 



3.3.1 Neutron-rich nuclei in the IsOd-lpOf shells 



Studies of extremely neutron-rich nuclei have revealed a number of intriguing new phenom- 
ena. Two sets of these nuclei that have received particular attention are those with neutron 
number N in the vicinity of the IsOd and O/7/2 shell closures (N ~ 20 and N « 28). Ex- 
perimental studies of neutron-rich Mg and Na isotopes indicate the onset of deformation, 
as well as the modification of the N = 20 shell gap for 32 Mg and nearby nuclei [j43 ]. In- 
spired by the rich set of phenomena occurring near the N = 20 shell closure when N ^> Z, 
attention has been directed to nuclei near the N = 28 (sub) shell closure for a number of S 
and Ar isotopes [0,0 where similar, but less dramatic, effects have been seen as well. In 
parallel with the experimental efforts, there have been several theoretical studies seeking to 
understand and, in some cases, predict properties of these unstable nuclei. Both mean-field 
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Table 1 

The computed and measured values of B{E2) for the nuclei in this study using e p = 1.5 and 
e„ = 0.5. 



B(E2; 0+ 2+) Expt B(E2, total) SM MC B(E2; 0+ 2 



22 Mg 

30 Ne 
32 Mg 

36 Ar 

40 g 

42 S 
42 Si 

44g 
44 Ti 

46 Ar 



458 ± 183 



454 ± 78 @3| 
296.56 ± 28.3 J57| 



334 ± 36 |4| 
397 ± 63 [El 



314 ±88 |5| 



610 ± 150 58 



196 ±39 p| 



334 


± 


27 


303 


± 


32 


494 


± 


44 


174 


± 


48 


270 


± 


66 


194 


± 


64 


445 


± 


62 


274 


± 


68 


692 


± 


63 


369 


± 


77 



342 |gl,171 |H 



448 |g],205 |H 



398 ||],390 |5]]] 
372 (431,465 |5 



260 p] 



271 g|,390 |5| 
460 |1,455 (5| 



46| , |47| and shell-model calculations p^j45|j48| -|52| have been proposed. The latter require 



a severe truncation to achieve tractable model spaces, since the successful description of 
these nuclei involves active nucleons in both the IsOd- and the lpO /-shells. The natural 
basis for the problem is therefore the full IsOd-lpOf space, which puts it out of reach of 
exact diagonalization on current hardware! . 



Shell-Model Monte Carlo (SMMC) methods [?, |53|j54| offer an alternative to direct diag- 
onalization when the bases become very large. Though SMMC provides limited detailed 
spectroscopic information, it can predict, with good accuracy, overall nuclear properties 
such as masses, total strengths, strength distributions, and deformation, precisely those 
quantities probed by the recent experiments. 

There is limited experimental information about the highly unstable, neutron-rich nuclei 
under consideration. In many cases only the mass, excitation energy of the first excited 
state, the B(E2) to that state, and the /3-decay rate is known, and not even all of this 
information is available in some cases. From the measured B(E2), an estimate of the 
nuclear deformation parameter, fa, has been obtained via the usual relation 



fa = 4n^B(E2; 0+ -> 2+)/3Zi? 2 e (78) 
with R = 1.2 A 1 / 3 fm and B(E2) given in e 2 fm 4 . 

Much of the interest in the region stems from the unexpectedly large values of the deduced 



For a treatment of the cm. problem, see Ref. [42| 
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02) results which suggest the onset of deformation and have led to speculations about 
the vanishing of the N = 20 and N = 28 shell gaps. The lowering in energy of the 
2f state supports this interpretation. The most thoroughly studied case, and the one 
which most convincingly demonstrates these phenomena, is 32 Mg with its extremely large 
B(E2) = 454 ± 78e 2 fm 4 and corresponding /3 2 — 0.513 |]4"3"|| ; however, a word of caution is 
necessary when deciding on the basis of this limited information that we are in the presence 
of well-deformed rotors: for 22 Mg, we would obtain fa — 0.67, even more spectacular, and 
for 12 C, j3 2 = 0.8, well above the superdeformed bands. 



Most of the measured observables can be calculated within the SMMC framework. It is 
well known that in deformed nuclei the total B(E2) strength is almost saturated by the 
0+ — > 2\ transition (typically 80% to 90% of the strength lies in this transition). Thus the 
total strength calculated by SMMC should only slightly overestimate the strength of the 
measured transition. In Table 1 the SMMC computed values of B(E2, total) are compared 
both to the experimental B(E2;0g S — > 2\) values and to the values found in various 
truncated shell-model calculations. Reasonable agreement with experimental data across 
the space is obtained when one chooses effective charges of e p = 1.5 and e n = 0.5. All of the 
theoretical calculations require excitations to the lpO /-shell before reasonable values can 
be obtained. We note a general agreement among all calculations of the B(E2) for 46 Ar, 
although they are typically larger than experimental data would suggest. We also note a 
somewhat lower value of the B(E2) in this calculation as compared to experiment and 
other theoretical calculations in the case of 42 S. Table 2 gives selected occupation numbers 
for the nuclei considered. We first note a difficulty in extrapolating some of the occupations 
where the number of particles is nearly zero. This leads to a systematic error bar that we 
estimate at ±0.2 for all occupations shown, while the statistical error bar is quoted in the 
table. The extrapolations for occupation numbers were principally linear. Table 2 shows 
that 22 Me; remains as an almost pure s<i-shell nucleus, as expected. We also see that the 
protons in 30 Ne, 32 Mg, and 42 Si are almost entirely confined to the sd shell. This latter is 
a pleasing result in at least two regards. First, it shows that the interaction does not mix 
the two shells to an unrealistically large extent. Second, if spurious cm. contamination 
were a severe problem, we would expect to see a larger proton O/7/2 population for these 
nuclei due to the Oc/5/2-0/7/2 "transition" mediated by the center-of-mass creation operator. 
The fact that there is little proton / 7 / 2 occupation for these nuclei confirms that the cm. 
contamination is under reasonable control. See Ref. [4~2~| for further details. An interesting 
feature of Table 2 lies in the neutron occupations of the N = 20 nuclei ( 30 Ne and 32 Mg) and 
the N = 28 nuclei ( 42 Si, 44 S, and 46 Ar). The neutron occupations of the two N = 20 nuclei 
are quite similar, confirming the finding of Fukunishi et al. and Poves and Retamosa 
f49fl that the iV = 20 shell gap is modified. In fact, the neutron O/7/2 orbital contains 
approximately two particles before the N = 20 closure, thus behaving like an intruder 
single-particle state. Furthermore, we see that 2p-2h excitations dominate although higher 
excitations also play some role. We also see that the neutrons occupying the lp0/-shell 
in N = 20 systems are principally confined to the O/7/2 sub-shell. The conclusions that 
follow from looking at nuclei with N > 20, particularly those with N = 28, are that the 
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Table 2 

The calculated SMMC neutron and proton occupation numbers for the sd shell, the 0/ 7 / 2 sub- 
shell, and the remaining orbitals of the pf shell. The statistical errors are given for linear extrap- 
olations. A systematic error of ±0.2 should also be included. The first row represents neutron 
results, while the second row represents protons. 





N,Z 


IsOd 


OA/2 


lpO/g/2 


22 Mg 


10,12 


3.93 ± 0.02 


0.1 ±0.02 


-0.05 ±0.01 






2.04 ±0.02 


0.00 ± 0.01 


-0.05 ±0.01 


30 Ne 


20,10 


9.95 ± 0.03 


2.32 ± 0.03 


-0.26 ±0.02 






2.03 ± 0.02 


-0.01 ±0.01 


-0.02 ±0.01 


32 Mg 


20,12 


9.84 ± 0.03 


2.37 ±0.03 


-0.21 ±0.02 






3.99 ± 0.03 


0.05 ± 0.02 


-0.05 ±0.01 


36 Ar 


18,18 


9.07 ±0.03 


1.08 ± 0.02 


-0.15 ±0.02 






9.07 ±0.03 


1.08 ±0.02 


-0.15 ±0.02 


40 g 


24,16 


11.00 ±0.03 


5.00 ± 0.03 


-0.01 ±0.02 






7.57 ±0.04 


0.54 ±0.02 


-0.12 ±0.02 


42 Si 


28,14 


11.77 ±0.02 


7.34 ± 0.02 


0.90 ± 0.03 






5.79 ± 0.03 


0.25 ±0.02 


-0.07 ±0.01 


42 S 


26,16 


11.41 ±0.02 


6.33 ±0.02 


0.25 ±0.03 






7.49 ± 0.03 


0.58 ± 0.02 


-0.09 ± 0.02 


44g 


28,16 


11.74 ±0.02 


7.18 ±0.02 


1.06 ±0.03 






7.54 ± 0.03 


0.56 ± 0.02 


-0.12 ±0.02 


44 Xi 


22,22 


10.42 ± 0.03 


3.58 ±0.02 


0.00 ± 0.02 






10.42 ± 0.03 


3.58 ± 0.02 


0.00 ± 0.02 


46 Ar 


28,18 


11.64 ±0.02 


7.13 ±0.02 


1.23 ±0.03 






8.74 ± 0.03 


1.34 ±0.02 


-0.08 ± 0.02 



A^ = 20 shell is nearly completely closed at this point, and that the A^ = 28 closure shell is 
reasonably robust, although approximately one neutron occupies the upper part of the lpOf 
shell. Coupling of the protons with the low-lying neutron excitations probably accounts for 
the relatively large B(E2), without the need of invoking rotational behavior. In Table 3 
we show the SMMC total Gamow- Teller (GT~) strength. We compare our results to those 
of previous truncated calculations, where available. In all cases, our results are slightly 
smaller than, but in good accord with, other calculations. Since we do not calculate the 
strength function, we do not compute /3-decay lifetimes. 
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Table 3 

The calculated total Gamow- Teller strength, GT~ , from this study. The results of other studies, 
when available, are presented for comparison. 



Nucleus 


SMMC 


Other 


22 Me 


578 ± 06 




30 Ne 


29 41 ± 25 




32 Mrr 


24 00 ± 34 




36 Ar 


2 13 ± 61 




40 g 


22.19 ± 0.44 


22.87 [51 


42g 


28.13 ±0.42 


28.89 [51 


42 Si 


40.61 ± 0.34 




44g 


34.59 ± 0.39 


34.93 [51 


44 Ti 


4.64 ± 0.66 




46 Ar 


29.07 ±0.44 


28.84 [51 



3.3.2 Heavy Sn isotopes 

Effective two-hole matrix elements are calculated based onaZ = 50, N = 82 asymmetric 
core and with the active P-space for holes based on the 2sy 2 , ld 5 / 2 , ld 3 / 2 , 0g 7 / 2 and 0hn/ 2 
hole orbits. The CD-Bonn model for the NN interaction was employed 0] and all diagrams 
through third-order in the G-matrix were included in the evaluation of the Q-box. Folded 
diagrams were again included through the method exposed above, see Refs. p||61| for further 



details. The corresponding single-hole energies £(^3/2) = 0.00 MeV, e(h 



els 



1/2J 



3/2) 

0.332 MeV, e(d+ /2 ) = 1.655 MeV and e{g, 



ni/2) 



0.242 MeV, 



7/2 j 



2.434 MeV are taken from 



Ref. [§2J and the shell model calculation amounts to studying valence neutron holes outside 
this core. The shell model problem requires the solution of a real symmetric n x n matrix 
eigenvalue equation 



H\9 k ) = E h \9 k ). 



(79) 



where for the present cases the dimension of the P-space reaches n w 2 x 10 7 . At present 
our basic approach in finding solutions to Eq. ( |T9"D is the Lanczos algorithm; an iterative 
method which gives the solution of the lowest eigenstates. This method was already applied 
to nuclear physics problems by Whitehead et al. in 1977. The technique is described in 
detail in Ref. [[S3]] , see also Ref. [Q. The results of the shell model calculation are presented 
in Table 4. All experimental information in the present analysis is taken from the data base 
of the National Nuclear Data Center at Brookhaven |53] . 

The isotopes above 116 Sn (heavy Sn) are treated based on the asymmetric Z = 50, iV = 82 
core. This simplifies the shell model calculation, but in addition it is of interest to see how 
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Table 4 

Exitation spectra for the heavy Sn isotopes. 



130 Sn 

r Exp. J 71 " Theory 


128 Sn 

J 71 " Exp. r Theory 


(2+) 1.22 2+ 1.46 
(4+) 2.00 4+ 2.39 
(6+) 2.26 6+ 2.64 


(2+) 1.17 2+ 1.28 
(4+) 2.00 4+ 2.18 
(6+) 2.38 6+ 2.53 


126 Sn 

r Exp. J w Theory 


124 Sn 

r Exp. r Theory 


2+ 1.14 2+ 1.21 
4+ 2.05 4+ 2.21 
6+ 2.61 


2+ 1.13 2+ 1.17 
4+ 2.10 4+ 2.26 
6+ 2.70 


122 Sn 

r Exp. J 71 " Theory 


120 Sn 

r Exp. J 71 " Theory 


2+ 1.14 2+ 1.15 
4+ 2.14 4+ 2.30 
6+ 2.56 6+ 2.78 


2+ 1.17 2+ 1.14 
4+ 2.19 4+ 2.30 
6+ 2.86 


118 Sn 

r Exp. r Theory 


116 Sn 

r Exp. J 71 " Theory 


2+ 1.22 2+ 1.15 


2+ 1.30 2+ 1.17 



successful a hole-hole effective interaction calculated with respect to Sn is. 

Only some selected states are displayed. First of all, the well-known near constant + — 
2 + spacing is well reproduced, all the way down to 116 Sn. Also the additional calculated 
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states are in very good agreement with experiment. However more detailed analysis of the 
results close to 116 Sn indicates that our effective two-particle interaction has difficulties in 
reproducing the shell closure which is believed to occur in this region. The increase of the 
the + — 2 + splitting is not as sharp as found experimentally, even if the phenomenon is 
rather weak in the case of Sn. We have observed a similar feature around 48 Ca Q which 
is generally agreed to be a good closed shell nucleus. There the deviation between theory 
and experiment is severe. Preliminary analysis indicates that our effective interaction may 
be slightly too actractive when the two particles occupy different single-particle orbits. 
This may be related to the radial wave functions which in our calculation are chosen to be 
harmonic oscillator functions. 



3.4 Inclusion of hole-hole contributions and single-particle propagators 



With the G-matrix defined according to the double-partitioned scheme we can easily solve 
Eq. (|27D through matrix inversion. The number of hole-hole and particle-particle configu- 
rations is then rather small, typically smaller than ~ 100, and a matrix inversion is then 
rather trivial. Before we discuss the solution of Eq. (|27|), it is always instructive to consider 
the contributions to second order in perturbation theory, i.e., diagrams (a) and (b) of Fig. 
[|. The external legs can be particle states or hole states. Diagram (a) reads 

(a) = \ E V SIj 7—T-: V v3^ (80) 

Z pq S bp bq 



and 

(h)-ly* v [12] 1 

\°) — 9 2^ V 12a/3J 1 , V a/3MJ- 

1 a g —S + E a + Bp 



We note here the minus sign in the energy denominator, since in the latter expression we 
are using the hole-hole term of the propagator of Eq. (^). If we use a double-partitioned 
G-matrix for say 16 and are interested in an effective valence space interaction for the IsOd- 
shell @ , then typically the single-particle orbits of the intermediate states will be represented 
by states in the lpOf major shell. Hole states are then defined by single-particle states in the 
0s and Op shells. Clearly, the number of two-body intermediate states is rather limited. To 
third order we have diagrams like (c) and (d) of Fig. |j. Diagram (c) is just the third-order 
equivalent of Eq. (|80D and reads 



(c ) - I V V [12] 1 1 ydi2] (82) 

pqrw " 



This means that the labels 1234 will refer to particle states in the lsO<i-shell. 
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while diagram (d) contains both a two-particle and a two-hole intermediate state and reads 



- 7 Yj ^Sj — — ; — V pqa(3j — v iphj- ( 83 ) 



Thus, solving Eq. (|27| ) will then yield contributions to the effective interaction such as the 
above expressions. 

Here we have also tacitly assumed that the energy denominators do not diverge, i.e., we 
have chosen an energy s so that we avoid the poles. This has always been the standard 
approach in calculations of shell-model effective interactions. To give an example, consider 
now diagram (b) and suppose that we are using harmonic oscillator wave functions. Let us 
also assume that the two hole states are from the Op-shell and that the valence particles are 
in the lsOrf-shell. If we rescale the energies of the valence space to zero, then the two-hole 
state would yield —28 MeV with an oscillator parameter b = 1.72 fm. If s = —28, the 
denominator diverges. In this case it is rather easy to obtain the imaginary part, and even 
if we were to chose s different from —28 MeV, the imaginary part will influence the real 



part of the effective interaction through dispersion relations, see e.g., Refs. |27||28|j68|] . It 
is therefore at best just a first approximation to neglect the imaginary term. Moreover, if 
we solve Dyson's equation for the self-energy, the single-particle energies may contain an 
imaginary part. Technically it is however not difficult to deal with imaginary contributions, 
one needs to invert a complex matrix rather than a real one. However, care must be 



exercised in localizing poles, see e.g., Ref. p6[ for a computational approach to this problem. 



In our actual calculations we will also follow Ref. . These technicalities will however be 



described elsewhere 57 . 



Using the double-partitioned G-matrix, we can then rewrite Eq. ( p7|) as 

r [ 12 l ( S )-G [12] +lV(7 [121 g [12] T [12] (s ) ( M ) 

1 1234jl 6 J — ^1234,7 + r, ^1256^ 1 5634,A 6 J; 

Z 56 

where G' 12 ^ is just the double-partitioned G-matrix discussed above. It is also energy de- 
pendent, in contrast to V. In case we were to employ this equation for effective interactions 
in the IsOfi-shell, the intermediate two-particle states would then come from just e.g., the 
lpO/-shell. This equation, which now is solved within a much smaller space than the origi- 
nal one spanned by the total Q pp , allows clearly for computationally amenable solutions. It 
corresponds to the so-called model-space approach to the solution of the Feynman-Galitskii 
equations as advocated by Kuo and co-workers, see e.g., Ref. for more details. Thus, 
a possible approach would consist of the following steps 



(i) Solve the G- matrix equation from Eq. (IS) using the double-partitioning scheme 



(ii) The next step is then to solve Eq. fl27|) and Dyson's equation for the self-energy. 

(iii) This scheme is iterated till self-consistency is achieved, see the discussion below. 
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We will however not employ this model-space scheme in our actual calculations. There are 
several reasons for not doing so. 



Let us first assume that we omit the [13] and [14] channels in our iterative scheme for Eq. 
(53). The next iteration of Eq. fl34|) would then look like 

p[12] _ p[12] [12] ^ [12] p [12] 

1 (1) - 1 (0) + 1 (0) y 1 (1) ' l 8b J 

[1 21 /rrr—n 

where the vertex function F, \ is the solution of Eq. (|8"4]). However, we cannot define the 
"bare" vertex T^3 to be the solution of Eq. fl34|) simply because then we would be double- 

M 2] _____ 

counting contributions. Thus, Tr Q \ has to equal the G-matrix. The only change in Eq. 
arises from the solution of Dyson's equation and thereby new single-particle energies. 

Let us then for the sake of simplicity assume that the single-particle energies are just the 
Hartree-Fock solutions. The problem we are aiming at arises at the Hartree-Fock level. 
In order to obtain Hartree-Fock solutions which are independent of the chosen harmonic 
oscillator parameter b, we typically need to include single-particle orbits from quite many 
major shells. Typical constraints we have found when we do so-called Brueckner-Hartree- 
Fock (BHF) calculations for finite nuclei is that we need at least 2n + I < 20 in order to 
obtain a result which is independent of the chosen fei value. The way we solve the BHF 
equations is to expand the new single-particle wave functions ip\, with A representing the 
quantum numbers nlj, in terms of harmonic oscillator wave functions, i.e., 

2n+Z<20 

IV*>= E~ <^I<U (86) 

a=l 

where cj) a are the harmonic oscillator wave functions with quantum numbers a = nlj and 
C are the coefficients to be varied in the Hartree-Fock calculations. The single-particle 
energies at the Hartree-Fock level are just 

e Q = ^« + E G(£ a + \ath) , (87) 

h 

where the single-particle states are just those of the harmonic oscillator. The (7-matrix used 
in the first iteration in the BHF calculation is the one given by the solution of Eq. fl48|) . The 
coefficients C a can then be obtained by diagonalizing a matrix of dimension N x N, where 
iV is the number of single-particle orbits with the same Ij values. As an example, suppose 
that we are considering the S\/2 state with / = and j = 1/2. With the above requirement 
2n + 1 < 20, we may have that A^ max = 10, the dimensionality being equal to the quantum 
number n. The way to proceed in a BHF calculation is to calculate the reference G-matrix 



8 Throughout this work our unperturbed single-particle basis par excellence will always be that 
of the harmonic oscillator. 
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Gf in Eq. ( f49|) once and for alia . Thereafter, the change in single-particle wave functions 
is introduced in the calculations of Eq. 



We see then that if we choose to do the Hartree-Fock self-consistency employing the double- 



partitioned G-matrix and summing the pphh diagrams as in Eq. (85), our single-particle 
basis will just be defined by the Os, Op, IsOd and IpOf shells in the case of an effec- 
tive interaction in the lsOd-shell. This is simply not enough in order to obtain a stable 
Brueckner-Hartree-Fock result. 

The reader could infer that why do we not perform first a BHF calculation for the G- 
matrix in Eq. (0), and then solve Eq. ( pT|) and stop there. This would be in line with 



the abovementioned model-space approach of Kuo and co-workers [gfj. However, if one 
calculates the self-energy by only including particle-particle intermediate states in the 
vertex function T, which is the case if we do a standard BHF calculation, one may seriously 
violate various sum rules, as demonstrated in Refs. p8| , |69| . Thus, to respect sum rules such 
as the conservation of number of particles constrains severely the way we solve Eqs. ( p7|) 
and (|25|). The vertex function needs both particle-particle and hole- hole intermediate states 
in order to eventually satisfy e.g., the conservation of the number of particles. 

Moreover, our interest lies in solving the Parquet equations. This entails simply that we 
perform the above self-consistency. It is therefore not only a matter of many-body aestethics 
for embarking on the solution of the Parquet equations. We are also able to satisfy various 
sum rules. In section [| we will also show that at every level of approximation, the solution 
of the equations in the [13] and [14] channels will also result in an antisymmetric vertex 
function in these channels. 

In section [5] we will come back to the technical solution of Eq. ( PTf ) and its iterations. 



4 Screening corrections and vertex renormalization, the equations for the [13] 
and [14] channels 



We start as in the previous section with the definition of the interaction vertices in the [13] 
and [14] channels and the corresponding integral equations. Thereafter, we discuss various 
approximations to these equations such as the summation of TDA and RPA diagrams. 
Eventually, the aim is to merge the discussion in this section and the preceeding one into 
equations for a self-consistent scheme which combines all three channels, namely the so- 
called set of Parquet equations to be discussed in section 



9 This matrix is typically set up in the relative and cm. system and calculated only once and 
for all for 6 = 1 fm. Other b- values involve simply a multiplication with a constant. The matrix 
Gf can therefore be used for other mass areas as well. 
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The equations for the renormalized vertex in the [13] and [14] channels have the same form 
as Eq. fl26D , namely 



[131 



^[13] + y[13] ( ^ )F [13] ) 



and 



r [14] = y[14] + y[14]^ F [H]_ 



(88) 



i9) 



The matrix elements which enter are however defined differently and the irreducible dia- 
grams of and can obviously not be the same. With irreducible in the [13] channel 
we will mean a diagram, which by cutting an internal particle-hole pair, cannot be sep- 
arated into a piece containing the external legs 1,3 and another piece containing 2,4 as 
external legs. The definition for the irreducible vertex in the [14] channel is similar and we 
illustrate these differences in Fig. Diagram (a) is just the lowest-order interaction in the 



(a) 



(b) 



(c) 



(d) 



Fig. 9. Examples of irreducible and reducible diagrams in the [13] and [14] channels. See text for 
further details. 

[13] channel and is therefore irreducible. Diagram (b) is an irreducible diagram in the [13] 
channel, whereas it is reducible in the [14] channel. Diagram (c) is in turn irreducible in 
the [14] channel and reducible in the [13] channel. Diagram (d) is an example of a diagram 
which is irreducible in both channels. This diagram stems from the [12] channel. 



The energy variables in these channels are, following Fig. || and Eqs. (15) and (|16|). 

t = e 3 - £i = e 2 - e 4 , (90) 
for the [13] channel and 

u = £i - e A = e 3 - e 2 , (91) 
for the [14] channel. Defining the unperturbed particle-hole propagators in the energy 



Qph 



representation as |36 

£[13] 

and 

<J[141 



t - e p + e h + it] t + e p - e h - it] 



(92) 



u - E, p + e h + i7] u + e p — e h — it] 



(93) 
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we arrive at the following equations for the interaction vertex in these two channels 

r£L(*) = vSSu + E vMJMr%j(t), (94) 

ph 



and 



T^LAu) = vSSu - E V^J^r%j(u). (95) 

ph 



These equations, together with Eq. (|27|) , can then form the basis for the first iteration in 
a self-consistent scheme for renormalization corrections of the Parquet type. The origin 
of the minus sign in Eq. (|95"D follows from the diagram rules |34j] and will be examplified 



below. A graphical view of these equations is given in Fig. The reader should also keep 
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Fig. 10. (a) shows the structure of the integral equation for the interaction vertex in the [13] 
channel, (b) represents the integral channel for the [14] channel. The coupling order is displayed 
as well. 



in mind the two contributions to the particle propagators of Eqs. 
36| for a physical interpretation. 



and f|93D. See Ref. 



In this section we will omit a discussion of the self-energy corrections which arise from 
these channels. This will be relegated to the next section. In the nuclear case, due to 
the strongly repulsive short-range character of the interaction, we will have to replace in 
actual calculations the bare interaction in the irreducible vertices [13] and [14] with the 
G-matrix discussed in the previous section. Doing this entails already a first step towards 
the Parquet set of equations, in the sense that we are including short-range correlations 
from the [12] channel. Here we will however limit the discussion to expressions in terms of 
the interaction V. The aim is to try to recover from Eqs. (0) and (|95|) the familiar TDA 
and RPA equations and the so-called self-consistent coupled equations of Kirson p5[. The 
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hope is that these intermediate steps can bridge the gap between the familiar TDA, RPA 
and G-matrix equations and the Parquet set of equations in section [5|. 



4-1 Screened ph and 2p2h interactions 



Here we study the screening of the particle-hole and the 2p2h interactions given in Fig. [I], 
indicated by V p h and V^h, respectively. Before we list the final expression, it is however 
instructive to consider the corrections to second order in the interaction V to the ph and 
2p2h vertices. In Fig. [TT| we display the second-order corrections to the ph diagrams of Fig. 




(a) 



(b) 



(c) 



(d) 



Fig. 11. Second-order perturbation theory corrections to the ph interaction vertex. 

|I|. Diagram (a) is the core-polarization correction term to the particle-hole interaction, and 
corresponds to a contribution from the [14] channel, as indicated by the coupling order. 
The term labeled (b) corresponds to the exchange term of (a) and is coupled in the [13]- 
order, see also the discussion in connection with Eqs. fl2~l~|) and ([22|). The other corrections, 
like (c) and (d) include particle-particle and hole-hole intermediate states, respectively. 
They are irreducible in both the [13] channel and the [14] channel, and can therefore enter 
the irreducible vertices of these two channels in later iterations. They are however not 
generated by various iterations of Eqs. fl94]) and fl95|). In fact, if we replace V by G in Eqs. 
(|9~4]) and (pop, diagram (c) is already accounted for by the G-matrix. It may however be 
included if the double-partitioned G-matrix of the previous section is used. Let us now 
look at the analytical expressions in an angular momentum coupled basis for diagrams 
(a) and (b) of Fig. [□]. Here we just include the first term of the propagators of Eqs. ( |9~2"D 
and (|93|). The second terms will give rise to the 2p2h contributions discussed below. In 
the following discussion we will also assume that the interaction V does not depend on 
the energy, although it is rather easy to generalize to an energy dependent interaction. 
Diagram (a) reads 



= -£(- 



U + S 



-v. 



[1-1] 



(96) 



7 



The factor (— ) 2j7 stems from the opening up and recoupling of an internal particle-hole 
pair |34| and the phase (— )>+J-r- J is needed in order to rewrite the matrix elements in 
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the coupling order of Eq. ([15]). The general structure of Eq. (|96| ) is just of the form 



—Vpl^Qph/ e ^Vph\ with e' 14 ' = e q + £ 7 — ep — e r = u + £ 7 — e r and we have defined 



u 



for the on-shell energy case. This is the equivalent of the energy variable of Eq. 
[12] channel. Diagram (b) is in turn given by 



(97) 



in the 



(b) = E(-) JV+J ^ J (-) 2J7 < 3 L™^ 



[13] 



r-y 



[13] /3r<ryJ' 



(9£ 



and we note that the contributions are clearly different. The minus sign in Eq. (|96| ) stems 
from the standard diagram rules [34j. In our use of the diagram rules below, we will omit 
the use of the rule for the number of external valence hole lines. In our case then, as can 
also be deduced from inspection of Fig. |ll|, diagram (a) has zero closed loops and three 
hole lines, giving thereby rise to a minus sign. Diagram (b) has an additional closed loop 
and thereby yielding the plus sign. The energy denominator is in this case 



J13] 



t + £- 



with 



(99) 



(100) 



We notice, using the relations discussed in Eqs. fl2~T|) and (|2~2"|), that diagram (a) is simply the 
exchange diagram of (b). We need however to include both diagrams in order to obtain an 
antisymmetric equation for the particle-hole channels which exhibits the same properties as 
the [12] channel shown in Eq. fl24|) . This is actually crucial in solving the Parquet equations. 
We wish namely that every iteration, with a given approximation to the vertex function V , 
preserves the antisymmetry property. This point cannot be emphasized enough. Let us 
now see what happens to third order in the interaction. Third order corrections to the ph 
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5^ 





V 




(a) 



(b) 



Fig. 12. Corrections beyond second order in the interaction V to the ph interaction vertex, (a) is 
in the [14] channel and (b) is in [13] channel. 



vertices (a) and (b) involving only ph intermediate states are shown in (a) and (b) of Fig. 
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|12| , respectively. The analytical expression for the third-order contribution (a) is given by 



1 ' ' ' " X] / VpyJaj — V}ss\j ~ . ~ _ I K/3q5 J 5 (10!) 



rs7(5 u + £ 7 - £ r 1 u + e s -e s 

with / = (— yr+j s +j-,+jg-2j ^_yj~f+2j s _ rp^jg e q Ua ti on has the general structure 

ir [14] Vph IA [14] Vph v [14] 
"ph e [ 14 ] v ph e [ 14 ] v ph • 

A similar expression applies to diagram (b), whose expression is 

(b) = E /<l/7— l -—V%— <L- (102) 

rs7l5 t + £ 7 - e r t + e<5 - £ s 

It has the general structure 

n [13] n [13] 
y[13] Vph ^[13] Vph_ v [13] 
Ph £ [13] P h e [13] P h ' 

But these expressions have the same sign! Diagram (a) counts now 4 hole lines, and (b) 
counts also 4 hole lines and 2 closed loops. However, there are three interaction terms V, 
and taking the exchange term of each of these in diagram (a) leads to the desired results, 
namely (a) = —(h), as it should. Thus, to third order we keep the antisymmetry property 
of T in the [13] and [14] channels. It is easy to see that in the [14] channel we will always 
have an alternating sign in front of each contribution, since every new order in perturbation 
theory brings a new hole line and no closed loop, and thus a new minus sign. In the [13] 
channel we have always one new hole line and one new closed loop for every new vertex. 
If we consider only the screening of the ph vertex, we can then set up a perturbative 
expansion in terms of the ph vertex for the vertex functions I^ 13 ' and I^ 14 '. For notational 
economy we will skip the Pauli operators Q^j | hp in the discussions below. It will always 
be understood that the intermediate states are two-body particle-hole states, |ph) or |hp). 
Consider e.g., 

r [14] _ y[U] _ y[U] J_y[14] , yM J_yM _L V ^ _ + ( W3 ) 
1 - "ph "ph e [i4] "ph + "ph e [i4] "ph e [i4] "ph "I" • • • , l-LUJj 

which can be summed up to yield 

p[14] _ y[U] _ y[U] 1 y[U] _ y[U] _ y[U] J_p[14] (104:) 

1 - "ph "ph [M] [ M ] V -V "ph e [ M ] i ' l iU4 J 

^ "ph 

which is the standard TDA expression for the ph term. The corresponding expression in 
the [13] channel results in 

^ 3] -v^ + v^r^. (105) 
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The signs agree with the expressions of Blaizot and Ripka p8| , see chapter 15 and Eq. 
(15.50). The summations in both channels ensures that the final vertex is antisymmetric 
and the combination of the latter two equations results in the familiar TDA equations, see 
e.g., Ref. [|36| for a matrix equation version. We next look at the 2p2h matrix element and 



show the corresponding corrections to second order in perturbation theory in Fig. 13. If 



(a) 



(b) 





(d) (e) (f) 

Fig. 13. Corrections to second order in V of the 2p2h vertex. 

we omit diagrams (e) and (f) which contain 2p and 2h intermediate states generated by 
the solutions in [12] channel, we have for diagram (a) 



£(-) >+jW (-) 2 ^! 



i 



■V. 



[14] 



(106) 



with the general structure 



-V, 



[14] 



V, 



[14] 



2p2h £ [i4] v 2p2h- 



(107) 



Note well the minus sign in front of u. The contribution from the propagator can in this 
case be retraced to the second term of the propagator of Eq. (|93|). Diagram (a) follows the 
coupling order of the [14] channel. It is also easy to see that diagram (b) is given by 



w = -E(-) jV+j ^ J (-) 2j7 d 



•v. 



14] 



and has the structure 

*ph £ [14] K 2p2h- 

Similarly, if we now move to the [13] channel we have the following expressions 



(106 



(109) 



(c) = £(- 

fry 



J * £ 



(3"/qrJ 



—t + e. 



■V, 



[13] 



(110) 
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and 



(d) = EH jV+jW (-)^CLj _, ■ P — V S$tJ> ( 1U ) 
with the general structure 

J_y[13] / 112 N 

K 2p2h e [i 3 ] * / 2p2h' l ±±z J 

and 

^[13] J_y [13] / 113 x 

P h e [13] V 2p2h; l i±0 J 

respectively. Diagram (c) is just the exchange of (a) and includes two 2p2h vertices, while 
diagram (d) is the exchange of diagram (b) and includes a ph vertex multiplied with a 2p2h 
vertex. We note again that the antisymmetry is ensured at a given order in the interaction 
only if we include the corrections at the same level in both channels. 

One can then easily sum up higher-order corrections to the 2p2h diagrams as well in both 
channels. The inclusion of the backward going particle-hole pair in the propagators of Eqs. 
(|92|) and (j93f ) ensures thus that we will also sum to infinite order 2p2h corrections. This 
leads ultimately to the familiar RPA equations, see e.g., Refs. p6|j3l 



A closer inspection of Eqs. ( |108| ) and fllUT ) shows that if we only include ph vertices, we 
could resum these corrections to infinite order for the 2p2h vertex by observing that the 
structure of such diagrams would be of the form (e.g., in the [14] channel ) 

p[14] _y[U] _ V [U]J_y[M] V [U]J_ V [M]J_ V [U] ( , 

1 2p2h — v 2p2h v ph [14] v 2p2h ^ v ph [14] v ph [14] v 2p2lT T • • • , \ LL ^J 



which can be summed up to yield 

1 

V 



p[14] _y[U) y[U) 1 y[U) ( -.-.^ 

1 2p2h _ v 2p2h v ph ri/|1 [141 2p2h' \ LLO ) 



and similarly for the [13] channel, but with a plus sign. The modification discussed in Eqs. 
( |104j) and ( |115|) serve to modify the propagation of a particle-hole pair and have normally 
been termed for propagator renormalizations, as can easily be seen from Eqs. (|104|) and 



( |115|) where the propagation of a free particle-hole pair is modified by the presence of the 
interaction V in the energy denominator. Other important processes which can affect e.g., 
various polarization terms are those represented by so-called vertex renormalizations, a 
term originally introduced by Kirson and Zamick [ff5| |. These authors studied the renor- 



malizations of the 2plh and 2hlp vertices as well, see also Kirson [25] and Ellis and Osnes 



26| for further discussions. We will therefore end the discussion in this section by looking 



at such renormalizations. 
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4-2 Further renormalizations 



In the previous subsection we dealt mainly with what has conventionally been labelled for 
propagator renormalizations. We will therefore extend the standard TDA and RPA scheme 
by looking at further ways of renormalizing interaction vertices. The approach discussed 
here follows Kirson 



25 1 . Extensions were made later by Ellis and Goodin JTTJ] and Ellis 



Mavrommatis and Miither [ 3~2| . We will limit the discussion here to the scheme of Kirson. 
We start therefore with the contributions to second order to the 2plh vertexS. These 




(b) 




(c) 





Fig. 14. The corrections to second order in V of the 2plh vertex. 



contributions are shown in Fig. |14[ Diagram (a) consists of a 2plh vertex multiplied with a 
ph vertex whereas (b) stems from the multiplication of a 2p2h vertex with a 2plh vertex. 
They both contain a particle-hole pair as an intermediate state and follow the coupling 
order of the [14] channel. The exchange diagram of (b) is given by (c), while that of (a) 
ig diagram (d). Diagrams (e) and (f) represent contributions from the [12] channel and 
are irreducible in both particle-hole channels. Before we sketch the general structure of 
the renormalization procedure of Kirson, it is instructive to consider again the equations 
to second order in perturbation theory, as the general expressions can be deduced from 
inspection of these diagrams. Diagram (a) is 



E( 

r7 



\jr+jj~J I 



)2i 7 y[14] 
I ' fyyraj 



U + E~ 



-v. 



[11] 



rsq"/Ji 



116) 



and the plus sign stems from the diagram rules ||34|| , i.e., we have two hole lines and no 
closed loop. The propagator is that arising from the first term in Eq. (|93|). The general 
structure is 



V, 



[14] 



V, 



[1-1] 



Ph ,[14] ^2plh- 



117) 
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The discussion here applies to the other interaction vertices discussed in Fig. |l|. 
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Diagram (b) reads 



(b) = E( 



\jr+j-y-J ( _\2j-yy[ 14 ] 



-U + £- 



'1181 



with the following structure 



[ii] 



[14] 



2plh e [14] ^2p2h' 

In this case the propagator stems from the second term in Eq. 
for e.g., the 2hlp vertices of Fig. |I[ 



(119) 

Similar equations arise 



Before we write down the self-consistent equations of Kirson |2q| , let us assume that we 
can approximate the 2plh vertex in the [14] channel by the first order term and diagrams 
(a) and (b). This renormalized vertex, which we here label V2 P ih, is then given by 



~ V [lA] + T/M—T/I 14 ] 4. V [lA] V [14] 

1p\h ~ v 2plh ' ph [14] v 2plh ' v 2plh [14] v 2p2h' 



(120) 



If we now allow for the screening to infinite order of the ph vertex given by Eq. ( |104|) 
and replace the 2p2h vertex in the above equation with Eq. ( |115| ) we obtain the following 
renormalization of the 2plh vertex 



T> _ T/[ 14 l i T/[ 14 l 



1 v m . v ^ 

,[14] _ T/ [14] 2plh 2plh 
^ph 



1 



[14] 



e [i4] _ V M V ^ 



'12V, 



A similar equation applies to the 2hlp vertex of Fig. [I] and for the [13] channel. Eqs. ( |104| ), 
( |115| ) and ( |121| ) form then the starting point for the approach of Kirson [^J. Examples 
of diagrams which can be obtained through the iterative solution of Eqs. (|104j) , (|115|) and 
(p.21[) are given in Fig. ITS. 





(a) (b) 

Fig. 15. Examples of diagrams which can arise from Kirson's self-consistent equations. 



The question now is however how to relate Eqs. (|104p , ( |115| ) and ( |121| ) with those from 
Eqs. ([94]) and (|95|) . This is rather trivial if we recall that the labels 1234 can, as was also 
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discussed in section 0, represent whatever single-particle states, either holes or particles. 
Thus, V1234 can represent a 2plh, 2hlp, 2p2h, 2p, 2h or a ph vertex. This means that, 



due to the choice of propagators in Eqs. fl92"|) and ( p3[ ) equations like Eq. ( |121| ) are already 
inherent in Eqs. (|94]) and d95|) . If we e.g., approximate Eq. (|)5|) to second order in the 
interaction V and let the single-particle labels 1234 represent a 2plh interaction vertex, we 
immediately reobtain Eq. (|120|) . If we let 1234 represent a 2p2h vertex, we find to second 
order diagrams (a)-(d) of Fig. ITS. 



Till now we have however refrained from discussing the contributions from the [12] chan- 
nel, examples were only shown in (c) and (d) of Fig. [11], (e) and (f ) of Figs. [13] and [14]. 
These diagrams cannot be generated by simply iterating the equations for the [13] and [14] 
channels, but they could enter as contributions to the irreducible vertex in the [13] and [14] 
channels from the first iteration in the [12]-channel. We see thus the emerging contour of 
an iterative scheme. The crucial point is however how to perform the next iteration of say 
Eqs. ( TO ), (1TT5D and ( pi]) or Eqs. 0, © and (|7|) from the [12] channel The question is 



how do we include the results from the first iteration into the next one, i.e., how to modify 
the bare vertices V^ 13 ' and in e.g., Eqs. (0) and (^) in order to obtain an effective 
interaction for the shell model. We have also not addressed how to deal with the solution 
of Dyson's equation for the one-body Green's function. We mention also that Ellis and 
Goodin |7l| included pp correlations, i.e., terms from the [12] channel such as diagram (e) 



of Fig. [TJ|, as well when they considered the screening of the 2plh and 2hlp vertices. Fur- 
thermore, as already mentioned in the introduction the authors of Refs. |32| , |33| extended 
the pp RPA to include the particle-hole (ph) RPA, though screening of the 2plh and 2hlp 



vertices was not included. In Ref. [Q however, a study with self-consistent single-particle 



energies was also performed. These works represent thus a first serious step towards the 
solution of the Parquet equations, i.e., a many-body scheme which solves self-consistently 
the equations in the [12], [13] and [14] channels, with the addition of the self-consistent 
evaluation of the self-energy. It ought also to be mentioned that one of the really first ap- 
plications for nuclear systems was performed in a series of papers by Dickhoff and Miither 



and co-workers |72[ for nuclear matter. These authors actually performed the first iteration 



of the three channels. 

Such a synthesis of the equations for the three channels discussed here will be made in the 
next section. 



5 Effective interactions for finite nuclei from Parquet diagrams 

The equations we discussed in the two previous sections can be generalized in matrix form 

as 

r = r m + r [ij] g [ij] r, (122) 
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where obviously [ij] represents a given channel, (/™ is the particle-particle, or hole- hole 
or particle-hole propagator. The propagator is a product of two single-particle propagators 
g which we do not specify any further here. They are defined by the solution of Dyson's 
equation in Eq. (P9"|). The irreducible vertices must appear in the solution of the self-energy, 
and conversely, the self-energy must appear in all single-particle propagators within the 
expressions for the three channels [12], [13] and [14]. 

For all of our practical purposes, the irreducible vertex used in all channels is the so-called 
free G-matrix defined in Eq. (|49|). We will explain more why this is our preferred choice 
below. Let us now define the contribution from the [12], [13] and [14] channels following 
Ref. [ I7[ . Eq. ([27]) is then rewritten as 

L = _ g f , (123) 

where obviously L stands for ladder and if we neglect hole-hole terms in Eq. (|2"T| ) we obtain 

L = AG = G - G F , (124) 

that is the ladder corrections beyond first order in our irreducible vertex function Gf, as 
discussed in section |[ The ladder term can then be rewritten as 

L = G F G [12] G F + G F G [12] L. (125) 

In a similar way we can define the diagrams beyond first order in the particle-hole channel 

as 

R [13] = p[13] _ (126) 

and 

R M = r [i4] _ Gjrj ( 127 ) 

where Gp now is coupled in either the [13] or [14] way. Rewriting where R refers to 
ring diagrams, we obtain 

rM = G F g [ij] G F + G F g [ij] R [ij] , (128) 

where ij stands for either [13] or [14]. The equation for the vertex function T in Eq. ( |122| ) 
becomes then 

r = G F + L + R^ +R^. (129) 

Note here that the vertex V can be represented in the coupling order of any of the above 
channels. Our convention is that of the [12] channel. 
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If Eqs. ( |125| ) and ( |128|) define the first iteration, it should be fairly obvious to see that the 
next iteration would be 

L = (G F + + §M (G F + + + (g f + Rim + R M) gl^L, (130) 

where the contributions i?' 1 ^ + -R' 14 ^ are recoupled according to the coupling order of the 
12-channel. These contributions are irreducible in the [12] channel. Similarly, for the rings 
we have 

rW = (p F + l + r^) (g f + l + rw) + [g f + l + r^) S [13] i? [ul , (131) 

and 

= {g f + l + r^) (g f + l + rM) + (g f + l + rM) g^r^. (132) 

Our scheme for calculating V will be an iterative one based on Eqs. (|129| )- (|132| ) and the 
solution of Dyson's equation for the single-particle propagator. This set of equations will 
then yield the two-body Parquet diagrams. Relating the above equations to the discussions 
of sections § and f|, it is rather easy to see that the G-matrix, TDA, RPA and Kirson's 
screening scheme are contained in Eqs. 



5.1 Petit Parquet 



The aim here is to present a numerically viable approach to the Parquet equations. Here 
we will limit ourself to just sketch the structure of the solution, more technical details will 



be represented elsewhere |67 



The iterative scheme starts with the solution of Eq. (|27|) . The bare vertex which is irre- 
ducible in all three channels is the so-called free G- matrix G F defined in Eq. (f49|) . Let us 
use the identity from Ref. |37j and rewrite the vertex function r' 12 ! as 



r [ 12 l =G F + G F 
where 



s - H s - H Q 



r [12] 



(133) 



G f 



V + V- 



-Gi 



We define then = Q pp + Qhh and use Eq. pTj ) to rewrite Q pp and obtain 



,[12] 
(o) 



G F — Gi 



Q 



hh 



Hr 



+ 



1 



-P^ 



^P- 



s-H P{s-H Q )- l P s-H 



r 



[12] 
(0) • 



(134) 
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The subscript (0) means that this is just the first iteration. The single-particle energies are 
the unperturbed harmonic oscillator energies. 



Since we will always deal with real single-particle energies, the two-body propagators of 
Eqs. (|28f) , fl92|) and fl93|) will have the same expressions except for the fact that the single- 
particle energies get renormalized after each iteration. This means that the denominators 
in Gfi rLy and subsequent iterations and the expressions for the ring diagrams from the 



[13] and [14] channels, can be rewritten via the simple relation 
1 1 

P- T m5{s). (135) 



s ±ir] s 

With this caveat, we can in turn obtain the real and imaginary parts of all matrices 
involved. 



The first step in our calculations is to evaluate Gp. It is calculated in momentum space for 
a series of starting energies (typically ~ 30). The principal value integral from Eq. ( |135|) 
is solved using Kowalski's [|73) method. This method ensures a numerical stable treatment 
with large numbers of mesh points in momentum space. From the principal value integral 
we can in turn define the real and imaginary part of Gp, see e.g., Ref. |28| for technical 
details. The matrix Gp is solved only once and transformed to a harmonic oscillator basis 
in the lab from the relative and center of mass system, see Ref. |4| for details. 

The second step is to solve Eq. ( |134| ), which is now a complex equation. Since the number of 
two-hole states is rather limited, typically < 100 for a given J even for heavy nuclei like Pb, 
the major problem in the matrix inversion of Eq. ( |134j ) resides in the dimensionality of P 
discussed in section [| What conditions our truncation of P and the number n 3 , is dictated 
by the Brueckner-Hartree-Fock (BHF) independence on the chosen oscillator parameter b. 
This choice leads to the inclusion of more than ten major shells in the computation of the 
G-matrix. For a given J value, the total number of two-body states needed can then be of 
the order ~ 10 4 . 

Our first approximation is therefore to truncate the available space of single-particle states 
to be within ~ 10 — 20 major shells. The choice being conditioned by the BHF criteria and 
the possibility to store these large matrices in the RAM of available computing facilities. 
If we choose to use a double-partitioned model space, we need to sum up further diagrams 
with particle-particle intermediate states. 

These two steps lead then to the first iteration of the ladders, i.e., 

L ( o) = rg 1 - Gp. (136) 

It contains both hole-hole and particle-particle intermediate states and is a complex matrix. 
The external single-particle legs can be particles or holes. Only unperturbed single-particle 
energies enter the definition of the two-body propagators. 
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The third step is to calculate the first iteration for the rings, namely 

R [ $ = (G F + L (0) ) G [13] (G f + L (0) ) + (G F + L {0) ) G [13] R [ $, (137) 

and 

R [ $ = (G F + L (0) ) g [14] (G F + L( 0) ) + (G F + L (0) ) G [li] R^. (138) 

The equations for L and R are all defined within a truncated Hilbert space. They can 
therefore be recast into matrix equations of finite dimensionality. Recall also that we need 
to recouple the contribution from the [12] into the relevant ones for the [13] and [14] 
channels. This is done employing Eqs. fll5|) and (|i~6|) . With these contributions, we can 
now obtain the vertex function T after the first interaction 

T(o) = G F + £(o) + -R(o) 1 + #(o) ■ (139) 

The fourth step is to compute the self-energy and thereby obtain new single-particle en- 
ergies. In so doing, care has to be exercised in order to avoid double-counting problems. 



A thourough discussion of this topic can be found in Ref. [H]. More details will also be 



presented in Ref. [ 67j ■ Through dispersion relations [28 we can in turn obtain the real part 



of the self-energy and our single-particle energies will be approximated by 

e a = t a + ReS Q . (140) 



This is our second approximation. The new single-particle wave functions of Eq. (p6|) are 
obtained by diagonalizing a matrix of dimension the quantum number n of 

the single-particle state a. 

The fifth step is to repeat steps 1-4 with the new single-particle energies till a predetermined 
self-consistency is obtained. But now the rings have to be included in all equations, i.e., 
we solve Eqs. dl29D-(p2|). 

The final vertex T can then be used to define a new effective interaction to be applied in 
shell model studies, where many more diagrams are considered than in present state of the 
art calculations, see e.g., Fig. 8 of Ref. |L7| for a list of diagrams to sixth order entering 
the definition of the irreducible vertex V. 



6 Perspectives 



The reader should always keep in mind that the many-body scheme we have focussed on 
is only one of several possible ways of calculating effective interactions. Other methods are 
also discussed in this volume. 
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The aim of this work has however been to show how one can practically implement the 
Parquet equations in order to obtain effective interactions for the nuclear shell model. The 
emphasis here has been to connect these equations with state of the art approaches to 
effective interactions. Thus, how to recover the standard G-matrix equation and the TDA 
and RPA equations from the Parquet equations has been outlined. Applications of this 
many-body scheme will be presented elsewhere HS7 |, although in depth discussions and 
applications of both the G-matrix, folded diagrams and other perturbative resummations 
can be found in e.g., Refs. @J2(|. There are also subtle technical details which deal with 
double-counting problems in the outlined iterative Parquet scheme that we have left out, 
due to space limits, in the discussion. We take the liberty here to refer to e.g., Refs. JETLIpJ 



They will also be discussed in Ref. 57 



In the scheme we have sketched, there are obviously other important many-body contribu- 
tions at the two-body level which cannot be generated by the two-body Parquet equations. 
We have also omitted any discussion of three-body terms. Such terms could be generated 
if we were to solve the three-body Parquet equations, see Lande and Smith [|17[]. An im- 
plementation of the three-body Parquet equations is in progress. Three-body terms, as 
also mentioned in the introduction, are believed to be important, it should just suffice 



to mention studies of the Triton [74]. With two-body forces only one is also not able to 



reproduce properly the nuclear matter saturation point |75| |. In connection with three 



body contributions one has to carefully distinguish between three-body forces and effective 
three-body interactions. To understand this point consider the following arguments from 
studies of infinite nuclear matter and finite nuclei. The last three years have seen quite 
some advances in many-body studies of dense infinite matter, see e.g., Refs. J75]-|77| for 
recent surveys. These results can be summarized as follows. Firstly, new NN interactions 
such as the CD-Bonn potential |5IJ, different Nijmegen interactions J7S| and the recent 
Argonne Vis interaction |79j all fit the set of scattering data of the Nijmegen group with a 
X 2 per datum close to 1. All these interactions, when applied to calculations of the equa- 
tion of state (EoS), yield essentially similar equations of state up to densities of 3 — 4 
nuclear matter saturation density for both pure neutron matter and /3-stable matter when 
the non-relativistic lowest-order Brueckner theory (LOB) is used. Other properties like the 
symmetry energy and proton fractions do also show a similar quantitative agreement, see 
Ref. |75[ for more details. Secondly, the inclusion of more complicated many-body terms 



at the two-body J75[ level does not alter this picture and even the recent summation of 
three-hole line diagrams of Baldo and co-workers |77| results in an EoS which is close to 
LOB when a continuous choice is used for the single-particle energies in matter [[f7]]. The 
latter are examples of effective three-body contributions. These findings are also in line with 
recent works on the energy of pure neutron drops, where three-body clusters are included 
30fl , and large-scale shell-model calculations of Sn isotopes including effective three-body 



interactions ||81|| . Differences do however occur when one introduces real three-body forces. 



These are necessary in order to reproduce the saturation properties of nuclear matter |75 
and the binding energy of light nuclei p| 
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The reader may then eventually ask why do we bother at all to solve the two-body set of 
Parquet equations when most likely these new effective interactions will not cure any of 
the problems seen in shell model studies or nuclear matter. One needs at least to include 
some realistic three-body force and thereby to solve the three-body Parquet equations. The 
problem however with real three-body forces is that presently we have no serious candi- 
dates which exhibit the same level of quality and sophistication as the NN interactions 
mentioned above. The introduction of such forces lead therefore to strong model dependen- 
cies. Moreover, to assess properly many-body terms at the two-body level is an important, 
and not yet solved, problem per se. The Parquet equations allow one also to include in a 
self-consistent way several many-body terms, fulfill crossing symmetries and satisfy certain 
sum rules. In addition, the Green's function formalism can be taylored for finite temper- 
atures. This is clearly of interest for studies of infinite matter at finite temperature, such 
as e.g., dense matter occuring in a newly born neutron star. The application of the Par- 
quet equation in nuclear and neutron matter will also allow for a consistent treatment of 
screening effects relevant for the pairing problem. Such work is in progress. 
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